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Abstract 

We present a study of diffractive associated J/tp + 7 production at the Fermilab 
Tevatron and LHC based on the Ingelman-Schlein model for hard diffractive scattering 
and the factorization formalism of NRQCD for quarkonia production, we find that 
this process {p + p — > p + J/V> + 7 + X) can be used to probe the gluon content of the 
Pomeron and test the assumption of diffractive hard scattering factorization. Using 
the renormalized Pomeron flux factor D ~ 0.11(0.052) , the single diffractive J/tp + 7 
production cross section at 4 < Pt < 10 GeV , — 1 < 77 < 1 region is found to be of 
the order of 3.0pb (8.5 pb ). The ratio of single diffractive to inclusive production is 
0.50%(0.15%) in central region at the Tevatron (LHC) for the gluon fraction in the 
Pomeron f g = 0.7, independent of the values of color-octet matrix elements. 



PACS number(s): 12.40Nn, 13.85.Ni, 14.40.Gx 



1 Introduction 



Early in the 1960's physicists already well realized that in high energy strong interaction the 
Regge trajectory with vacuum quantum number, the Pomeron, plays a particular and very 
important role in soft processes, such as the energy dependence of atot(s), the behavior of the 

i el 

elastic differential cross section ^^-at small the single and double diffractive dissociation 
processes in hadron-hadron collisions After the advent of QCD , physicists begin to 
study the nature of the pomeron in the QCD framework [[|. Since the Pomeron carries the 
quantum number of the vacuum, in QCD language, it is a colorless entity, which leads to the 
diffractive events are characterized by a large rapidity gap, a region in rapidity devoiding 
of hadronic energy flow, this distinct class of events are observed by the ZEUS and HI 
Collaboration at HERA in deep inelastic scattering region ||, which offer a unique chance 
to study the the soft scattering process with a hard virtual photon probe, therefore provide 
a more complete understanding of QCD . 

In an orginal paper |Q, It was suggested that hard diffractive scattering processes would 
give new and valuable insight about the nature of Pomeron. the Ingelman-Schlein (I-S) model 
assumes that the hard diffractive scattering processes can be calculated in a factorized way 
: first, a Pomeron is emitted from the diffractively scattered hadron, then partons of the 
Pomeron take part in hard subprocesses, the results will be diffractively produced high-Pp 
jets in P — P(P),7 — P collisions or larger rapidity gap events in diffractive deep inelastic 
scattering (DDIS) , so the partonic struction of the Pomeron could be established and stud- 
ied experimentally. These pridictions were confirmed subsequently. The UA8 Collaboration 
at CERN S ppS collider with y/s = 630 GeV studied the diffractive jet distribution , indi- 
cates a dominant hard partonic structure of the Pomeron, however, the diffractive dijet event 
topology alone can not distinguish between a hard-quark or a hard-gluon struction in the 
Pomeron [H. The DDIS and dijet photoproduction experiments at HERA have shed light on 
the partonic structure of the Pomeron. Combining the measurements of the diffractive struc- 
ture function in DDIS and the photoproduction jet cross sections, the ZEUS Collaboration 
gives the first experimental evidence for the gluon content of the Pomeron and determines 
that the hard gluon fraction of the pomeron f g is in the range 0.3 < f g < 0.8, independent 
of the validity of the momentum sum rule for the Pomeron and the normalization of the flux 
of Pomeron from the proton|| . The HI Collaboration also determines the fraction of the 
momentum of the Pomeron carried by the hard gluon, which is f g ~ 0.9 at Q 2 = 4.5 GeV 2 
and f g ~ 0.8 at Q 2 = 75 GeV 2 . The partonic structure of the Pomeron, based on the 
I-S model for hard diffraction , was also studied by the CDF Collaboration at the Tevatron 
through diffractive W production and dijet production |J, which give further evidence 
for the hard partonic stucture of the Pomeron. Combining these two experiments, the CDF 
Collaboration determined the hard-gluon content of the Pomeron to be f g = 0.7 ± 0.2, this 
result is also independent of the Pomeron flux normalization or the validity of the momentum 
sum rule for the Pomeron. 

However, there is a problem. The I-S model is based on the assumption of diffractive 
hard scattering factorization. Recently, a factorization theorem has been proven by Collins 
T0|| for the lepton induced diffractive hard scattering processes, such as DDIS and diffractive 
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direct photoproduction of jets. This factorization theorem justifies the analysis given by the 
ZEUS and HI Collaboration, and establishes the the universality of the diffractive parton 
distrubutions for those processes to which the theorem applies. In contrast, no factorization 
theorem has been established for diffractive hard scattering in hadron-hadron collisions. As 
already known before the advent of QCD, factorization fails for hard processes in diffractive 
hadron-hadron scattering 11]. Furthermore, in order to preserve the shapes of the M 2 and t 
distributions in soft single diffraction (SD) and predict correctly the experimentally observed 
SD cross section at all energy in PP collisions, Goulianos |12j has proposed to renormalize 



the Pomeron flux in an energy- dependent way, this approach indicates the breakdown of 
the triple- Regge theory for soft single diffractive exciation, and implies that diffractive hard 
factorization is likely breakdown in hadron-hadron collisions which appears to be confirmed 
by CDF experiments [|]. At large \t\ , where perturbative QCD applies to the Pomeron, it 
was proven that there is a leading twist contribution which broken the factorization theorem 
for diffractive hard scattering in hadron-hadron collision [T3). Evidence for this substantial 



coherent perturbative contribution has been observed by the UA8 experiment in diffractive 
jet production, in which the square of the proton's four-momentum transfer t is in the region 
—2 < t < — 1 GeV 2 ||. On the other hand, for the Pomeron at small \t\ , nonperturbative 
QCD phycics dominates, where the I-S model without much of a coherent contribution 
would be appropriate. Therefore, various hard diffractive processes may be considered to 
probe the partonic struction of the Pomeron, and test the hard diffractive factorization in 
hadron-hadron collisions |14] Recently, Alero et al. ||16|| extract the parton densities of 



the Pomeron from the HERA data on DDIS and diffractive photoproduction of jets, and use 
the fitted parton densities to predict the diffractive production of jets and weak bosons in 
PP collisions at the Tevatron. The predicted cross sections are substantially high than the 
experimental data, which signals a breakdown of hard scattering factorization in diffractive 
hadron-hadron collisions. 

In this paper, we will discuss another diffractive process, the associated J ftp + 7 single 
diffractive production at large Pt'- 

P + P^P + J/^ + j + X. (1) 

This process is of special significance because the produced large Pt J ftp is easy to detect 
through its leptonic decay modes and the J ftp's large Pt is balance by the associated high 
energy photon. We will see through the following calculations, although the production cross 
section of associated J/tp + j is sensitive to the color-octet matrix elements, the ratio of the 
single diffractive production cross section to the inclusive production cross section are not 
so. We find that the ratio is sensitive to the product Df g , where D is the renormalization 
factor of the Pomeron flux(which indicates the factorization broken effects, )and f g is the 
hard gluon fraction in the Pomeron. The measurement of this process at the Tevatron and 
LHC would shed light on the nature of the Pomeron and test the diffractive hard scattering 
factorization. Furthermore, this process is also interesting to the study of heavy quarkonium 
production mechanism. 

Our paper is organized as follows. We describe in detail our calculation scheme in Sec. II. 
A brief introduction of the heavy quarkonium production mechanism, the hard subprocesses 
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for associated J /if) + 7 production and a summary of the kinematics related to the J / if) Pt 
distribution are given in this section. Our results and discussions are given in Sec. III. 



2 Calculating Scheme 

Based on the I-S model for diffractive hard scattering, the associated single diffractive 

production process at large Pt consists of three steps(shown in Fig.l). First, the Pomeron is 
emitted from the proton with a small squared four- momentum transfer \t\ . Second, partons 
from the anti-proton and the Pomeron scatter in the hard subprocesses that produce the 
almost point-like large Pt pair of cc and a associated photon. In the third step, J /if) is 
produced from the point-like cc via soft gluon radiation. 



2.1 Heavy quarkonium production mechnism 

Prior to 1993, the conventional wisdom for heavy quarkonium production was based on the 
so-called color-singlet model (CSM) (17 ] which assumes that the heavy quark pair is produced 



in a color-singlet state with the right quantum numbers of the final heavy quarkonium in the 
hard subprocesses and the belief that the dominant production processes for color-singlet cc 
were the Feynman diagrams that were lowest order in a s . However, at the Tevatron, the 
CSM prediction for prompt J / if) production is an order of magnitude smaller than the CDF 
data at large Pt- The first major conceptual advance in heavy quarkonium production was 
the realization that fragmentation dominates at sufficiently large Pt [19] which indicates 



that most charmonium at large Pt is produced by the fragmentation of individual large Pt 
partons. The fragmentation functions are calculated in CSM , for example, the fragmen- 
tation fuction for g — > J /if) is calculated from the parton process g — > cc + gg in CSM . 
Including this fragmentation mechanism brings the theoretical predictions for prompt J/ if) 
production at Tevatron to within a factor of 3 of the data [p0[| . But the prediction for the 



if)' remains a factor of 30 below the data even after including the fragmentation contribution 
( the if) 1 "surplus" problem ). Furthermore, the presence of the logarithmic infrared diver- 
gences in the production cross sctions for P-wave charmonium states and the annihilation 
rate for Xcj — > qqg indicate CSM is incomplete. All these problems indicate that some 



important production mechanism beyond CSM needs to be included pl|[p^| . And so the 
color-octet mechanism is proposed which is based on the factorization formalism of nonrel- 
ativistic quantum chromodynamics (NRQCD) p4||[|25| . Contrary to the basic assumption of 



CSM, the heavy quark pair in a color-octet state can bind to form heavy quarkonium. 



Although the color-octet fragmentation picture of heavy quarkonium production |22| has 
provided valuable insight, the approximation that enter into fragmentation computations 
break down when a quarkonium's energy becomes comparable to its mass. The fragmentation 
predictions for charmonium production are therefore unrealiable at low Pt. In the case 
of T production, the exist data at the Tevatron are in the Pt < 15 GeV region, which 
significantly disagree with the fragmentation predictions. Based upon the above several 
recently developed ideas in heavy quarkonium physics, Cho and LeibovichPpT] identify a 
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large class of color-octet diagrams that mediate quarkonia production at all energies, which 
reduce to the dominant set of gluon fragmentation graphs in the high P? limit. By fitting the 
data of prompt J /if) and T production at the Tevatron, numerical values for the long distance 
matrix elements are extracted, which are generally consistent with NRQCD power scaling 



rules ||27|| . In order to convincingly establish the color-octet machanism, it is important to 
test whether the same matrix elements be able to explain heavy quarkonim production in 
other high energy processes, such as inclusive J/if> production in e + e~ annihilation on the 
Z° resonace and inelastic J/if) photoproduction at HERA et al. . 

2.2 The hard subprocesses for associated J/ifj + j production 

J I if) is described within the NRQCD framework in terms of Fock state decompositions as 

\J/i/j > = 0(1) |cc[ 3 5{ 1} ] > +0(v)\ccfpf ] ]g > + 

0(^)1^% > +0(v 2 )\cc[ 3 S?'%g > + 
0(v 2 )\cc[ 3 P^%g > + ■■-, (2) 

where the cc pairs are indicated within the square brackets in spectroscopic notation. The 
pairs' color states are indicate by singlet (1) or octet (8) superscripts. The color octet cc 
state can make a transition into a physical J /if) state by soft chromoelectric dipole (El) 
transition(s) or chromomagnitic Ml transition(s) 

(cc - )r+ l L (l,8) WW _ (3) 



NRQCD factorization scheme |2l has been established to systematically separate high and 



low energy scale interactions. It is based upon a bouble power series expasion in the strong 
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interaction fine structure constant a s = *p- and the small velocity parameter v. The pro- 
duction of a (cc)[ 25+1 Lj 1,8 ' ) ] pair with separation less than or of order — are perturba- 
tively computable. The long distance effects for the produced almost point-like cc to form 
the bound state are isolated into nonperturbative matrix elements. Furthermore, NRQCD 
power counting rules can be exploited to determine the dominant contributions to vari- 
ous quarkonium pro cesses [27]. For direct J /if) production, the color-octet matrix elements, 



< 0|O 8 J/V, [ 3 ^i]|0 >,<0\O J 8 /lp [ 1 S }p >,< 0|O^pPj]p > are all scaling as m 3 c v 7 c . So these 
color-octet contributions to J/ip production must be included for consistency. 

The partonic level subprocesses for associated J ' /ip + 7 production are composed of the 
gluon fusion subprocesses, which are sketched in Fig.2. These are 

9 + 9 - T + WM\'S®], 

g + g - 1 + (cc)[ 1 Si 8 \ 3 PP}. (4) 

The quark initiated subprocesses(gg channel) are strongly suppresed and will be neglected 
further. The color-singlet gluon-gluon fusion contribution to associated J /if) + 7 production 
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is well known [ESI: 



^(singlet) = ^ 



s 2 (s - Am 2 c ) 2 + t (t - 4m 2 ) 2 + u (u - Am 
(s - 4m 2 ) 2 (t - Aml) 2 (u - Ami) 2 



where 

s = (Pi + P2) 2 , i = (pa - P) 2 , « = (p x - P) 2 . 
The overall normalization M\ is defined as 
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where 



and e 2 = § 



Gi{J/il>) 



<o|o 1 w [ 3 5 1 ]|o>_r(j/^^ J +/- 



3m 2 . 



|-7re 2 a 2 



(5) 
(6) 
(7) 
(8) 



The average-squared amplitude of the subprocess g + g — > 7 + (cc)[ 3 5 , | 8 ' 1 ] can be obtained 
from the average-squared amplitude of g + g — > 7 + (cc)[ 3 S , j 1 ' ) ] by taking into account of 
defferent color factor. The result is 

da 



% ig + g->i + (cc)i 3 s[ 8) ]^ 1 + J/ij} 
at 

— 15 E|M( 5 + ^7 + (cc)[ 3 M 1) ])| 2 



<0|C>8^[ 3 5i]|0> 



167TS 2 6 

1 151 1 



24m r 



16vrs 2 6 8 12 



GA 4 2 2 2 

64^e c e m^, 



s 2 (s - 4m 2 ) 2 + t 2 (t - Am 2 c ) 2 + u 2 (u- Am 2 c ) 2 



<0\O^[ 3 S 1 }\0> . 



24m r 



(s - Aml) 2 (t - 4m 2 ) 2 (u - 4m 2 ) 2 



The average- squared amplitudes of the subprocesses g + g — ► 7 + (ccJ^Sq ] and g + g 
7 + (cc)[ 3 Pj 8 ^] can be found in p9 |, 



(9) 



at 



da 

di 



' E|Mfo + -> 7 + (cc)^])! 2 ■ J- <0|O 8 W [ 1 5 ]p >, 



167TS 2 ' 
[g + g^ 1 + {cc)fPf ) }^ 1 + J/^} 



^E^(S+S-T+ (cc)[ 3 Pf])| 2 • <0|O 8 ^[ 3 P ]p >, 



5m r 



where the heavy quark spin symmetry 

<0|O 8 W [ 3 P J ]p >= (2J + 1) <0|O 8 J/ ^[ 3 P ]p > 

is exploited. 



(10) 



(11) 
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2.3 The P T distribution of J/<if> 

Now we consider the Pt distribution of J/ip produced in process 

p{P p ) + p(P p ) - p(P' p ) + J/4>(P) + 7 (A0 + X (12) 

Based on the I-S model for diffractive hard scattering, the differential cross section can be 
expressed as 

da = fip/ p (£, t)/ s //p(xi, Q 2 )f g /p(x 2 , Q 2 ) ( ^d£dtdx 1 dx 2 di, (13) 

where £ is the momentum fraction of the proton carried by the Pomeron, t = (P p — P p ) 2 is 
the squared of the proton's four-momentum transfer. fip/ p (£,t) is the Pomeron flux factor 

d 2 asD/d£dt g(0) . 1 _ 2a(t) 2 , s 
= K^~ 2a ^F 2 (t), (14) 



where the parameters are choosen as JT2 



AT = 0.73Ge1/ 2 , a(t) = 1 + 0.115 + 0.26(GeV~ 2 )t, F 2 {t) = e 4M . (15) 

In the following calculation, we use the renormalized flux factor for the Pomeron, proposed 
by Goulianos JT3 in order to preserve the shapes of the M 2 and t distribution in soft single 
diffraction (SD) and predict correctly the observed SD cross section at all energies in pp 
collisions, 

fE!J p {Z,t) = Df J p /p (Z i t), (16) 
the renormalization factor D is defined as 

D = Min(l,^) (17) 

with 

N= dd dtf IP/p (U), (18) 

where £ m « n = ^ with Mq = 1.5GeV A2 (effective threshold) and £ max = 0.1(coherence limit). 
At the Tevatron energy (y/s = !800GeV), D = \. 

We now consider the kinematics. In the pp cm. frame, we can express the momenta of 
the incident p, p and gluons etc. as 

P P = ^(1, 1, 0) 

Pp = ^(i, - 1, 0) 
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Pip = ^(1, 1, 0) 
Pi = 1, 0) = x 1 P IP 



P2 = ^2(1, - 1, 0) = x 2 P p (19) 



where the first component is the energy, the second is longitudinal momentum, and the third 
is the transverse component of the four-momentum, x±,x 2 are the momentum fractions of 
the gluons. The momenta of the outgoing J/ifj are given by 

P=(E,P L ,P T ) = (E,P T shrj,P T ), (20) 

where Pt is the transverse momentum of J/ip, V is the pseudo-rapidity of J /if) and E = 

m \ + PT c h 2 V ■ 
The Mandelstam variables are defined as 

8 = (P p + Pp)\ 

s = (pi +P2) 2 = x x x 2 is, 
i = (pa -Pf =m\ - y/sx 2 (E + P T shri), 

u = ( Pl -P) 2 =ml-^x 1 ^(E-P T shr ] ). (21) 



Usings + t + u = mi, we have 



y/sfrijE - P T shrj) - m\ 
X2 x^s - ^~s{E - P T shrj) [ ' 



In order to obtain the distribution in the transverse momentum Pt for the process Eq.(|P2] 
we express the differential cross section as 



dcr = fp P l(U)f 9 /ip(xi,Q 2 )f g/P (x 2 ,Q 2 )^J(^^)d^dtdx l dr ] dP T , (23) 



where the Jacobian can obtain from Eqs.(^Tj) and Eq . ([22|) , 

j, x x x 2 i = 2sxiX 2 jP^chr] 
l xi^P T ~ E[xi£s - y/s(E + P T shrj)} ' 1 ' 

Then the Pt distribution of J /if) is expressed as 

^r= d v dd dxj dtf§,%(U)f 9/ ip(xuQ 2 )f g/ p(x 2 ,Q 2 )J(^)^, 

(25) 
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where the allowed regions of are given by 

E + P T sh V - 



Xi r 



£(V^ -E + P T sh V ) 



P + Prs/ir? - 

€dw = —j= . t-j — 7— • (26) 

V-s — P + PtsHt] 

In order to suppress the Reggon contributions, we set (, up = 0.05 as usual . 



3 Numerical results and discussions 

Now, we are ready to show the numberical results from the analytic expressions giving in 
the previous section. For numberical predictions we use m c = 1.5GeV, A4 = 235MeV, and 
set the factoriztion scale and the renormalization scale both equal to the transverse mass of 
J/ip, i.e., Q 2 = m 2 -. = (rr& + P|) . For the color-octet matrix elements <0|Cg^[ 3 S'i]|0>, < 

OlOg^pSyp > and <0|C 8 W [ 3 P ]P > we use the values determined by Beneke and Kramer 
| 3Cfl from fitting the direct J /if) production data at the Tevatron [31] using GRV LO parton 



distribution functions 



<0|O 8 J/,/ '[ 3 5 1 ]|0>= 1.12 x 10- 2 GeV 3 , 



<0|C 8 J/,/ '[ 1 5o]p > <0|C 8 W [ 3 P ]P >= 3.90 x 10~ 2 GeV 3 . (27) 



m 2 



Since the matrix elements < 0|(9 8 ^[ 1 S'o]p > and < 0|C 8 ^[ 3 Po]p > are not determined 
separately, we present the two extreme values allowed by Eq. (|27|) as 



l S saturated case : <0|C 8 /V '[ 1 ^o]P > = 3.90 x 10~ 2 GeV 3 

<0|O 8 W [ 3 P ]p> =0, 

3 Pj saturated case : <0\O^{ 3 P }fi > = 1.11 x lQ~ 2 mc 2 GeV 3 

<0|C 8 /V, [ 1 5 ]p > =0. (28) 

For the parton distribution functions, we use the GRV LO gluon distribution function for 



the anti-proton |32| and the hard gluon distribution fundtion for the Pomeron || 

xf 9 /ip(x, Q 2 ) = f g 6x(l -x), f g = 0.7 ± 0.2. (29) 

We use the central value of f g above for numberical calculation. We neglect any Q 2 evolution 
of the gluon density of the Pomeron at present stage. 

With all ingredients set as above , in Fig. 3 we show the Pp distribution Bj^ for single 
diffractive J/ip + 7 production in pp collisions at yfs = 1.8TeV , integrated over a pseudo- 
rapidity region — 1 < rj < 1 (central region). Where B = 0.0594 is the J/ip — > fi + n~ leptonic 
decay branching ratio. The lower solid line is the color-singlet gluon-gluon fusion contri- 
bution, the lower dashed line represents 1 Sq 8 ^ -saturated color-octet contribution, the lower 
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dotted one is 3 Pj 8 ' ) -saturated color-octet contribution and the lower dash-dot-dotted line is 
3 S[ 8 ^ color-octet contribution. For comparision, we also calculate the inclusive associated 
J /if} + 7 production p + p — > J/ip + '-f + X in the same kinematic region, the results are 
shown as the upper lines of Fig.3. The code for the lines are the same as the single diffractive 
production case. As shown in Fig.3 , the color-octet 3 Sf^ contribution is strongly suppressed 
compared with the others over the entire Pt region considered in both production cases. The 
1 Sq -saturated and 3 Pj 8 ^-saturated contribution are smaller than the singlet contribution 
where Pt < 5 GeV, though their contribution dominate in the high Pt region, their differen- 
tial cross section in the high Pt region is much small than that at low Pt region , so integrated 
over the Pt region, their contribution ( Bo = 0.08 pb and 16 pb for 1 S*Q 8 ^-saturated contribu- 
tion to SD and inclusive production respectively; Bo = 0.06 pb and 12 pb for 3 Pj 8 ' ) -saturated 
contribution to SD and inclusive production respectively ) are smaller than the color-singlet 
contributoin {Bo{singlet) = 0.10 pb and 20 pb for SD and inclusive production respectively 
). The total SD (inclusive) production cross section at y/s = 1.8 TeV, intergrated over 
4 < Pt < 10 Gev in central region is Bo SD {cen.) = 0.18 pb ( Bo Inclustve {cen^) = 36 pb ) for 
^-saturated case, Bo SD {cen.) = 0.16pb ( Bo Inclusive {cen.) = 32 pb ) for 3 Pj 8) -saturated 
case. The ratio of the total SD production cross section to that of inclusive production in 
the central region R cen - are 0.50% for both ^q 8 -* -saturated case and 3 Pj 8 ^-saturated case. 

In Fig. 4 , we show the Pt distribution of J /if) integrated over a pseudo-rapidity region— 4 < 
i] < — 2 (forward region), at y/s = 1.8TeV, the code of lines is the same as Fig.3. One can 
find the differential cross section is significantly smaller than that in central region. So 
the diffractively produced J /if} should be concentrated in the central detectors of collider. 



The same character was found early by Berger et al. [TJJ] in the rapidity distribution of 
SD production of bottom and top quarks. This contrasts with the naive expectation that 
diffractively producted system appera only at large rapidity. The total SD (inclusive) pro- 
duction cross section at y/s = 1.8 TeV, intergrated over 4 < P T < lOGev in forward 
region is Bo SD {fwd.) = 5.7 x 10~ 2 pb ( Bo Indusive {fwd.) = 23 pb ) for ^-saturated case, 
Bo SD {fwd.) = 5.1 x 10~ 2 pb ( Bo Inclusive {fwd.) = 20 pb ) for 3 Pj 8) -saturated case. The 
ratio of the total SD production cross section to that of inclusive production in the forward 
region are 0.25% for both ^-saturated case and 3 Pj 8 ' ) -saturated case. 

In TABLE I. , we show the ratio 

J SD J Inclusive 

R ^ - -k'^lPT- (30) 

at y/s = 1.8 TeV in central region and forward region, where 

do SD _ do SD {smglet) do SD { 3 S{ 8) ) dof^Sf) 
dP T ~ dP T + dP T + dP T ^ 

for 1 S , g 8 ' ) -saturated case, and 

do SD _ do SD {singlet) do SD { 3 S{ 8) ) do SD { 3 P^) 

775 1 775 ( 32 ) 



dP T dP T dP T dP 
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Table 1: The ratio R(Pr) as a function of Pt at the Tevatron in central region (—1 < f] < 
1) R cen -(P T ) and forward region (-4 < 77 < -2) R fwd \P T ). 



P T (GeV) 


4 


5 


6 


7 


8 


9 


10 


R cen (P T )(%) 


0.52 


0.52 


0.50 


0.48 


0.47 


0.46 


0.44 


Rl wd -{P T ){%) 


0.28 


0.28 


0.28 


0.28 


0.29 


0.29 


0.30 



for 3 Pj -saturated case. 

The ratio R(Pt) for 1 5 , Q 8 ' ) -saturated case is the same as that for 3 P® -saturated case. From 
this table, we can see that R(Pt) is almost constant for Pj- < 6 GeV both in central and 
forward regions. R(Pt) varies slowly as Pt increase. Because the differential cross section in 
the high Pt region is much small than that at low Pt region , integrated over the overall Pt 
region, the ratios R cen - and Rf wd - are not sensitive to the Pt smearing effects . Forthermore 
, we have varied the color-octet matrix elements <0\O^[ 3 Si]\0>,<0\O^[ 1 S Q \'Q > +f§ < 

0|Og^[ 3 Po]P > by multiplied them by a factor between ^ and 2, the result ratios R cen - 
and R fwd - is the same. This character demonstrates that the ratios R cen - and Rf wd - are 
not sensitive to the values of color-octet matrix elements. The ratios R cen - and Rf wd - are 
proportional to Df g , hence are sensitive to the gluon fraction of the Pomeron f g and the 
renormalization factor D which indicates the factorization broken effects. From other single 
diffractive production experiments, f g can be determined, as the CDF Collaboration at 
Tevotran and the ZEUS Collaboration at HERA DESY have done, the renormalization 
factor D can be determined precisely from those ratios, vice versa. So measuring these 
ratios can probe the gluon density in the Pomeron and shed light on the test of diffractive 
hard scattering factoriztion theorem. 

In Fig. 5, we show the Pt distribution of J/ip integrated over a pseudo-rapidity region— 1 < 
i] < 1 ( central ) at LHC energy y/s = 14TeV with the factor D = 0.052 calculated from 
Eq. (|17D . The code of lines is the same as Fig.3. As expected, the increasing of the differential 
cross section with cm energy is slowed down by the renormalization of the Pomeron flux. 
The total SD (inclusive) production cross section at a/s = 14 TeV, intergrated over 4 < Pt < 
lOGev in central region is Ba SD (cen.) = 0.51pb ( Ba Indusive (cen.) = 3.4 x 10 2 pb ) for 1 S$ 8) - 
saturated case, Ba SD (cen.) = 0.46 pb ( Ba Indusive (cen.) = 3.0 x 10 2 pb ) for 3 Pj 8) -saturated 
case. The ratio of the total SD production cross section to that of inclusive production in 
the central region R cen - are 0.15% for both -saturated case and 3 Pj 8 ^-saturated case. 

In the above calculations, we have set the values of factor D according to Eq.([T7|) , it is 
I at Tevatron energy and 0.052 at LHC energy. But the value is not unique, since it may 
change with different choice of the parameters such as M Q and £ ma x i n Eq. (|T8|) . If we use the 
central value of D = 0.18 ± 0.04 measured by the CDF Collaboration at the Tevatron, the 
predicted ratio in the forward region Rf wd - at the Tevatron , taking into account bothp, p can 
be diffractively scattered, will be 0.81% , which is close to the measured rate of diffractive 
dijet production with two jets of Et > 20GeV , 1.8 < \r)\ < 3.5 and 771772 > . 
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Experimentally, the nondiffractive background to the diffractive associated J /if) + 7 pro- 
duction must be dropped out in order to obtain useful information about the natrue of 
Pomeron and the factorization broken effects, this can be attained by performing the rapid 
gap enalysid ad in the CDF diffractive dijet experiment. 

In conclusion, in this paper we have shown that the diffractive associated J/if> + 7 pro- 
duction at large Pt is sensitive to the gluon content of the Pomeron and the factorization 
broken effects in hard diffractive scattering. Though the sigle diffractive and inclusive pro- 
duction cross section is sensitive to the values of coler-octet matrix elements, the ratio of 
single diffractive to inclusive J /if) + 7 production is not so and proportional to Df g , hence 
are sensitive to the gluon fraction of the Pomeron f g and the renormalization factor D which 
indicates the factorization broken effects. So experimental measurement of this ratio at the 
Tevatron and LHC can shed light on the nature of Pomeron and test the assumption of 
diffractive hard scattering fratorization. 
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Figure Captions 



Fig.l. Sketch diagram for diffractive associated + 7 production in the Ingelman-Schlein 
model for diffractive hard scattering. 

Fig. 2. The subprocess g± + g<i — > 7 + cc[ 2S+1 L^' 8 ^] — > 7 + J/^ in the NRQCD framework. 
Fig. 3. Transverse momentum of J/ip(Pr) distribution!? da/dPr , integrated over the J/ip 
pseudo-rapidity range \rj\ < 1 (central region), for single diffractive (lower) and inclusive 
(upper) associated J/ip + 7 production at the Tevatron. Here B is the branching ratio of 
J/ijj — > fi + fi~(B = 0.0594). The solid line is the color-singlet gluon-gluon fusion contri- 
bution, the dashed line represents 1 So 8 " 1 -saturated color-octet contribution, the dotted one 
is 3 Pj 8 ' ) -saturated color-octet contribution and the dot-dot-dashed line is 3 S[ 8 ^ color-octet 
contribution. 

Fig. 4. Transverse momentum of J/^(Pt) distribution!? da/dPr , integrated over the J/ip 
pseudo-rapidity range —4 < r\ < — 2 (forward region), for single diffractive (lower) and 
inclusive (upper) associated J/x/j + 7 production at the Tevatron. Here B is the branching 
ratio of J/ip — > fi + fi~(B = 0.0594). The solid line is the color-singlet gluon-gluon fusion 
contribution, the dashed line represents ^o^-saturated color-octet contribution, the dotted 
one is 3 Pj 8 ' ) -saturated color-octet contribution and the dot-dot-dashed line is 3 S[ 8 ^ color-octet 
contribution. 

Fig. 5. Transverse momentum of J/i(j(Pt) distribution!? rfa/rfPj' , integrated over the J/ip 
pseudo-rapidity range < 1 (central region), for single diffractive (lower) and inclu- 
sive (upper) associated J/ip + 7 production at LHC. Here B is the branching ratio of 
J/ip — > fi + fi~(B = 0.0594). The solid line is the color-singlet gluon-gluon fusion contri- 
bution, the dashed line represents 1 S , Q 8 ' ) -saturated color-octet contribution, the dotted one 
is 3 Pj 8 ' ) -saturated color-octet contribution and the dot-dot-dashed line is 3 S[ 8 ^ color-octet 
contribution. 
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